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ABSTRACT 

A new gauge-invariant approach for describing cosmological perturbations is developed. It is 
based on a physically motivated splitting of the stress-energy tensor of the perturbation into 
two parts — the bare perturbation and the complementary perturbation associated with stresses 
in the background gravitational field induced by the introduction of the bare perturbation. The 
complementary perturbation of the stress-energy tensor is explicitly singled out and taken to 
the left side of the perturbed Einstein equations so that the bare stress-energy tensor is the sole 
source for the perturbation of the metric tensor and both sides of these equations are gauge 
invariant with respect to infinitesimal coordinate transformations. For simplicity we analyze the 
perturbations of the spatially-flat Friedmann-Lemaitre-Robertson- Walker (FLRW) dust model. 
A cosmological gauge can be chosen such that the equations for the perturbations of the metric 
tensor are completely decoupled for the /locb^oi, and hij metric components and explicitly 
solvable in terms of retarded integrals. 
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1. Introduction 

The relativistic theory of perturbations in spatially homogeneous and isotropic cosmological models was first 
developed by E.M. Lifshitz (Lifshitz 1946; Lifshitz & Khalatnikov 1963). At present, two approaches are 
commonly used for the analysis of cosmological perturbations: the first approach is based on globally defined 
coordinates in the perturbed universe (Lifshitz 1946; Lifshitz & Khalatnikov 1963; Sachs & Wolfe 1967; Silk 
1968; Futamase 1989; Ma & Bertschinger 1995) and the second one (Bardeen 1980; Kodama & Sasaki 1984; 
Stewart 1990; Mukhanov, Feldman & Brandenberger 1992; Durrer 1990, 1994; Ma & Bertschinger 1995), 
generally known as the covariant approach, does not invoke any global coordinates but is based on a slicing 
and threading of the space-time with a specific choice of an orthonormal frame at each point. The main 
difficulty associated with the two approaches is that 8T a fj has always been used as the physical source for 
the primordial perturbation of the metric tensor h a p via the linearized Einstein equations. Such a direct 
interpretation of 8T a @ does not completely elucidate the physics underlying the perturbation theory. In 
general 5T a p must depend not only on the perturbation of matter but on the metric tensor perturbation h a $ 
as well. Nevertheless, the a priori dependence of 5T a p on h a p has remained implicit in previous approaches. 
It is the purpose of this paper to make this dependence explicit by developing a new approach to cosmological 
perturbation theory. 
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The theory developed here is based on a physically motivated splitting of the perturbed stress-energy 
tensor into two parts — the bare perturbation and the complementary perturbation associated with stresses in 
the background gravitational field induced by the introduction of the bare perturbation. The complementary 
perturbation of the stress-energy tensor is explicitly singled out and taken to the left side of the perturbed 
Einstein equations so that the bare stress-energy tensor turns out to be the sole source for the perturbation 
of the metric tensor. We also require that both sides of the linearized Einstein equations be independently 
invariant under gauge (Lie) transformations of the metric tensor perturbations. This approach allows us 
to find a new cosmological gauge such that the perturbation equations decouple and their solutions can be 
found explicitly in terms of retarded integrals. 

The development of the new approach was independently started in papers by de Vega, Ramirez & 
Sanchez (1999) and Kopeikin (2000), where a new cosmological gauge in the theory of cosmological pertur- 
bations was explored. Our aim here is to generalize the treatment started in these papers for the case of the 
dust-dominated spatially-flat FLRW cosmological model, which is adopted here for the sake of simplicity. 
Possible extension to other models with different equations of state of the background matter can be pursued 
in a similar way. 

We describe the background universe in Section 2. Then, in Section 3 the gauge invariance and in- 
tegrability of the Einstein equations are used for a unique scalar-tensor decomposition of the perturbed 
background stress-energy tensor in such a way that the Einstein equations for the metric perturbations are 
explicitly obtained with the bare stress-energy tensor of matter acting as the source of the perturbation of 
the metric tensor. A new cosmological gauge is discussed in Section 4 for reducing the gauge freedom of the 
linearized Einstein equations and bringing them to the form of the wave-type equations which are explicitly 
solved in Section 5 in the same way as d'Alembert's equation is solved in flat space-time. 

In what follows we use the 'geometrized' units in which G = c = 1 and the conventions adopted in the 
textbook of Misner, Thorne & Wheeler (1973). Greek indices in the full Einstein equations are raised and 
lowered with the help of the complete metric g a p and those in the linearized Einstein equations are raised 
and lowered with the help of the background metric g a[j . 



2. The background universe 

Let us choose the background cosmological model to be the spatially-flat FLRW space-time given by the 
metric 

9 a0 = a 2 ( V )f al3 , ^ = _L_ r /3, (1) 

where hereafter the overbar refers to the background, f a p = diag(— 1, +1, +1, +1) is the Minkowski metric, 
n is a dimensionless temporal coordinate related to the cosmic time t by the first-order ordinary differential 
equation dt = a(iq)dri, and a(r/) is a scale factor with the dimension of length. The cosmic time coincides 
with the proper time of static observers in the background space-time. The background Einstein equations 
read 

G a p = Ra/3 - -jg a pR = 87rT Q) 3 , (2) 

where R a p is the Ricci tensor, R = R° a is the scalar curvature, T a p = p u a U/3 is the stress-energy tensor 
of the background matter, and p(r]) and u a (n) are the density and four- velocity of the matter, respectively. 
The background universe is spatially homogeneous and isotropic, hence, u a {v) = a vl)Voa = ~ a ( r l)^a- Let 
us introduce the Hubble 'parameter' 7^(77) = a(rj) / a 2 (rj) , where the overdot denotes the time derivative with 
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respect to 77, so that 

R a p = ^ {g a p - 2u a u p ) + 3H 2 g ap , R = Q^L + 2H 2 ^j . (3) 

Equations (l)-(3) have a unique solution (Landau & Lifshitz 1971): 

a(v) = ^, H( V ) = W , P(V) = ^, t= — , (4) 

where 77 = 770 = 1 at the present epoch, and Ho is the present value of the Hubble parameter H(r]o). It is 
worth noting for the following calculations that H = — {2>/2)H 2 a and u a \p ~ HP a p, where P a p = u a up + g a p 
is the projection tensor on the hypersurface orthogonal to the unperturbed four- velocity u a . 



3. Basic assumptions 

Let us assume that the background space-time metric g af} is weakly perturbed by the presence of a dis- 
turbance with the stress-energy tensor Tfffl of arbitrary origin (e.g., a galaxy, a background fluid density 
perturbation, a cosmic string, etc.). This bare perturbation, in the absence of any interaction with the 
background matter, would be expected to move in the background space-time in such a way that 

T ( r )/3 „ = , (5) 

where the vertical bar denotes a covariant derivative with respect to the background metric g a/3 . The 
tensor must be understood as a weak perturbation of the background space-time metric only; however, 

equation (5) does not prohibit the energy density contrast /Too to be very large in some local region 
or even be singular. On the other hand, if one wants to consider the evolution of primordial cosmological 
perturbations that are not so well localized such as, e.g. a black hole, it must be assumed that the averaged 
energy density contrast {T^)/ (Too) <C 1 in most regions of the space-time manifold. 

The presence of the bare perturbation implies that the full space-time metric g a p can be written as a 
linear sum of the background metric g a/3 and the space-time perturbation h a p 

g a p{n,x) =g a p + h af }(ri,x) , (6) 

where we have introduced a four-dimensional coordinate chart x a — (77, x) on the background space-time 
manifold. In what follows it is convenient to introduce a new variable defined as 

VW3 = Kp - ^ g a[i h , (7) 

where h = h a a = g al3 h a p = —tp = —g at3 ip a p. The Einstein equations for the full metric read 

Ga/3 = Ra/3 ~ T^japR = 8vrT Q/ 3 . (8) 

These equations can be linearized after making use of equations (2) and (6). The result is 

G a p = G a p + SG a p , T a p = T a p + 5T a p . (9) 
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Here ST a p is a perturbation of the background stress-energy tensor and 5G a p is a perturbation of the Einstein 
tensor that in the case of homogeneous and isotropic background space-time can be written as 

SG afi = ~ (^ap ]v \ v + g a pB v ]v -B a \f}- B Pla j + m v {a i\) p)v - hi^ aP - 1 (li afi - ^g af3 ~Rj V , (10) 

where B a = ^>a \v Assuming that the background Einstein equations (2) are valid, the linearized Einstein 
equations read 

SG afj = 8tt ST a0 . (11) 

It is worth noting that the bare perturbation interacts gravitationally with the background space- 
time causing perturbations both of the background geometry and the background stress-energy tensor. As a 
result of this gravitational interaction a complementary stress-energy tensor is induced which explicitly 
depends on the metric perturbation h a p. For this reason in a linear approximation, the dressed perturbation 
of the background stress-energy tensor 8T a p can be written as a linear sum 

6T aP =Tif+T<$. (12) 

Because the complementary perturbation depends on h a p we can take to the left side of the linearized 
Einstein equations (11). The bare perturbation remains in the right side of equation (11) and has now a 
clear physical meaning: its introduction into the background manifold originates the deviation from the 
background metric g a p and leads to the perturbed space-time manifold. The complementary perturbation of 
the background stress-energy tensor is also generated by the bare perturbation and, hence, can be represented 
as 

T aj3 = u a u/3(Sp + Sp) + -pu a 5u fj + pu l3 5u a + g afj Sp , (13) 

where Sp, Sp, and Su a are as yet undetermined perturbations of the background energy density pressure, and 
the four- velocity respectively. It is worth noting that though the background pressure in the dust-dominated 
cosmological model is identically equal to zero, its perturbation 5p would be associated with the stresses of 
the background gravitational field and could be different from zero; that is, the bare perturbation disturbs 
the Hubble flow thereby inducing a purely gravitational effective pressure Sp. 

(c) 

The structure of these perturbations and, as a consequence, that of the complementary tensor T^J is 
completely determined by two main constraints. The first constraint comes from the integrability condition 
for equation (11) that is a direct consequence of the linearized Bianchi identity for the full Einstein equations 
(8). Accounting for equation (5) in linear approximation the Bianchi identity reads 

Ti c)u W + T JT^ - T JY v a0 = , (14) 

where Ta a p and ST^g are perturbations of the Christoffel symbols 

SKp = \(h» aW + h^ a -hJ v ) . (15) 

The second constraint on the structure of comes from the gauge invariance of the Einstein equations 
and the bare perturbation. That is, if one chooses a slightly different coordinate chart 

x' a =x a -C(ri,x), (16) 

the gauge (Lie) transformations 

SG' a0 {ri, x) = SG a0 (r,, x) + C^G af3 (v) , ST^( V , x) = 5T afi {n, x) + C^T a p( v ) , (17) 
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where denotes a Lie derivative along the vector field £, must preserve the form of the linearized Einstein 
equations (11) 

6G' a/3 (r,,x)=8ir8T^(r 1 ,x) . (18) 

Here 8G' aj3 and 5T' al3 are respectively defined by exactly the same relations as 5G a p and 8T a p after making 
the replacement h a p — ► h' a p — h a p + t, a \f} + £/3|a- We emphasize that in the linear approximation the bare 
stress-energy tensor remains invariant under the gauge transformations (16), i.e. T'^\r],x) = T^\r],x), 
since in our approach this tensor is the source for the perturbations of both the background geometry and 
the background stress-energy tensor. Hence, we conclude that in the new gauge T^(n, x) = T^(r), x) + 
CgT a p(i]). We can now use the results of the present section to determine the structure of the complementary 
stress-energy tensor. 



4. Gauge-invariant structure of the complementary stress-energy tensor 

The second constraint on the structure of the T^, discussed in Section 3, requires the introduction of sup- 
plementary fields $(77, x) and Z a (r), x). We assume in analogy with equation (6) for the metric perturbations 
that these scalar and vector fields can be expanded around their background values 

$=$(ri) + <l>(ri,x), Z a =( a (v) + ( a (r 1 ,x) . (19) 

Generalizing ideas developed in Lifshitz (1946) and (Bardeen 1980; Kodama & Sasaki 1984; Stewart 1990; 
Mukhanov, Feldman & Brandenberger 1992), we assume that the background matter perturbations can be 
linearly expressed in terms of the tensor field ip a 0(v> x )> the vector field ( a (r], x ), and the scalar field <j>{r], x) 
with coefficients that are proportional to the background quantities depending only on time n. We demand 
that under the Lie transformations (Weinberg 1972) 

i>' a p(v,x) = i> a p(v,x)+£ a \0 + £0\ a -9 a p^ ltl , (20) 

c a (v,x) = uv,x) + c al ,e + c^e la , (21) 

4>'(ri,x) = cb( V ,x)+^e, (22) 

our linear expressions for the matter perturbations in terms of ipa/3, Ccn an d 4> must change in accordance 
with the general gauge transformations 

Sp'^x) = Sp(ri,x)+p^e , (23) 

6p'(t],x) = Sp(r],x), (24) 

5u' a (r),x) = 5u a (r],x)+u a \^ + %£f a , (25) 

for our dust model. 

Straightforward calculations based on these consistency requirements reveal that one can choose the 
gauge vector field Z a = and <f> = — 3 In a(n), so that (j>\ a (r]) = 3Hu a . Furthermore, the matter perturbations 
can then be represented in the following scalar-tensor form 

6p = ^(V^-ItV) -^-g^^, (26) 
S P = UT^^-^TA+^n-^^, (27) 
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pSu a = ~ (T^W - 1 TV^ m q + £ 0, Q + 0|„) . (28) 

The complementary stress-energy tensor is now determined by equations (13) and (26)-(28), and is 
given as a sum of two pieces depending separately on the metric tensor perturbation ip a p and the scalar 
perturbation <fi, 

71(c) _ rpW , rp(4>) f Q\ 

T { $ = \g a ,{T^^ v -\Ti>j , (30) 

lif = \g a pT ^+^{u a ^ [3 +u f} 4>\ a -g ap u»4>^) . (31) 
The gauge transformation of the complementary tensor is given by 

T' a f( V , x) = T^(rj, x) + £^T a/3 (r?) , (32) 

where 

C^T^M = T apw T + T av C\p + Tp„ T| Q (33) 

is a Lie derivative of the background stress-energy tensor along the vector field £ Q that is generated by 
the coordinate transformation (16). Taking into account the gauge transformation equations (17) and (18), 
and the fact that in the linear theory the bare stress-energy tensor does not change under the gauge 

transformations (20)-(22), we conclude that the linearized Einstein equations can be written in an arbitrary 
coordinate system as 

8G a0 -8nT$ = 8nT$ , (34) 

where 5G a p is defined by equation (10). 

The supplementary scalar field <p is determined by making use of the first constraint on T^, i.e. the 
integrability condition (14), which gives a specific equation relating the vector functions B a and the scalar 
field </>; this equation is invariant with respect to the gauge transformations (20)-(22) and has the following 
form 

u a B a = -\e u a u^ afj - - \u<t> + ^4> W \ a ■ (35) 

We interpret this relation as the second-order differential equation for the scalar field <fi whose solutions 
depend on our choice of the (gauge) vector functions B a . 

The difference between our treatment of the perturbation of the stress-energy tensor 5T a p and the 
gauge-invariant approach adopted by previous authors (Bardeen 1980; Kodama & Sasaki 1984; Stewart 1990; 
Mukhanov, Fcldman & Brandcnberger 1992) is that we clearly distinguish the a priori bare perturbation from 
the complementary perturbation induced by the introduction of the bare perturbation into the background 
space-time. In addition, we are able to represent the complementary perturbation of the background 
stress-energy tensor as an explicit function of the tensor and scalar perturbations that obey well-defined 
hyperbolic differential equations. 

Our approach to cosmological perturbations appears to be rather suitable for situations where the bare 
perturbation can be thought of as an isolated physical system, since for such systems can be immediately 
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determined explicitly as, e.g., in the case of cosmic strings or domain walls (Kibble 1976; Vilenkin 1985). It 
would be interesting to explore in detail the relationship between our treatment and previous approaches. 
This is beyond the scope of the present endeavor, however, and will be dealt with in a future study. Here 
we only delineate some of the novel features of our formalism. 



5. The gauge-invariant linearized Einstein equations for gravitational perturbations 

Making use of formulas (29)-(31) for the complementary tensor T^, one can transform the linearized Einstein 
equations (34) to the form 

+ 9a0 BV \v -Ba\P - B t3\a -^R V (cc^l3) v + ^Ri>af}+ (R a ~ T^apRj 4> 

+ V a f3 (-R^Vv _ R 0) + 2H <t>\l3 + "/3 4>\ a - Vafl ^<?V) = -16ttT^ , (36) 

such that the left side is gauge invariant. This can be proved by applying the gauge transformations (20)— (22) 
to the left side of Einstein's equations (36) and accounting for the fact that the bare stress-energy tensor is 
gauge invariant by definition. Equation (35) for the scalar field <j) can be re-written in the source-free form 

<l>\ - \ h2 +\* + Su^a/j) - 3Itt*B a = . (37) 

This equation is also invariant under the gauge transformations (20)-(22). 

The set of equations (36)-(37) constitutes the basis for the new gauge- invariant approach to cosmological 
perturbations of the dust-dominated FLRW model. We have weighty arguments that an analogous approach 
can be worked out for the more general 'canonical' equation of state of the background matter p = ap with 
arbitrary numerical value of a as well as for the case of a universe with A-term and a background fluid with 
components that have different equations of state. 

The vector field B a and the scalar field <p are arbitrary gauge functions related through equation (37). 
A specific choice of these functions restricts the gauge freedom and can simplify the linearized Einstein 
equations. After choosing a specific gauge, the residual coordinate freedom is defined by the functions 
describing the gauge transformations (16). These functions obey the following inhomogeneous equation 

C\ + Tp e = B' a ( V , x) B a ( V , x) , (38) 

that is obtained by differentiation of equation (20) and making use of the definitions B' a = ip a \p and 

B a = i)^ a \p- Once B a is specified, four restrictions are imposed on the choice of the coordinate system. 
Then, equations (38) describe four residual degrees of gauge freedom in this coordinate system. Thus, 
the total number of functional restrictions on the eleven degrees of freedom of the gravitational field ip a fj 
and the scalar field <j) is eight, which means that we have three independent variables in the case of a free 
gravitational field propagating on the curved background instead of two that would be expected in the case 
of asymptotically-flat space-time. Two of the three degrees of freedom describe <g> and © polarizations of 
gravitational waves and the third degree of freedom belongs to the scalar field (f>. 
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6. New cosmological gauge 



Our new gauge-invariant approach to the theory of cosmological perturbations allows us to decouple the 
linearized Einstein equations. This simplifies the task of finding their solutions. Such a decoupling can be 
achieved in the framework of a new cosmological gauge that is defined by the condition 



(39) 



leading to the second-order hyperbolic equation for the scalar field 

J« ■ 3 



+ -H 2 

l« 2 



and a particular solution of equation (37) for the vector field B a , 

B a = 2H (u a 4> - • 
It is convenient to introduce the new metric variables (p a p'- 

ip a p(v, x) = a 2 (r]) tp a p(r), x) , <p(r},x) = f al3 ip afj (r],x) . (42) 
Then, Einstein's equations (36), in the gauge defined by equation (41), can be reduced to the following form 



(40) 



(41) 



U<p af > - 2H<p a0 , o +H 2 [6° a <p 0o + 6°^ a0 + 5%{ip - 20)] = -16rr T £ 



(m) 
a/3 



(43) 



where □ = f at3 d a dp, and H — a/a. Equations (40) and (43) are equivalent to the following set of inhomo- 
geneous hyperbolic equations 





_Adx 
77 dr] 
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7] dij 


rf ' 






4 cVoo 
■q dr] 






4 dipoi 
rj dr) 


4(^ i 
r] A 


= -l&rT^. 




n<p i:i - 


4 dipij 
r] dr] 


= -16^"\ 



_rp(m) 



8x 
rf 



(44) 
(45) 
(46) 
(47) 
(48) 



where \ = ^00 - 4> + </?/ 2 - 

The residual gauge freedom is restricted by the equations 



4 ^° , 4£° 
r] dr] 



+ 



if 



0, 



4<9f 

7/ i9r7 







(49) 



which can be derived from equation (38) after making use of equations (41), (20), and (22). It is evident that 
the new gauge (39) allows us to make time transformations along with transformations of spatial coordinates 
using solutions of equations (49). Thus, in the case of free gravitational waves one can construct an analog 
of the transverse-traceless (TT) gauge widely used for the description of physical effects of gravitational 
radiation in asympotically-flat space-times (see, e.g., Landau & Lifshitz (1971), Weinberg (1972), Misncr, 
Thorne & Wheeler (1973) as well as Grishchuk (2001) and references therein). 
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7. Solutions of the linearized Einstein equations 



All differential equations given in the previous section have the following symbolic form 



ad b i 

2 



-4775(77, a;) , 



(50) 



where a and & are constant numerical coefficients. This equation can be solved for arbitrary a and b using 
standard techniques (Arfkcn & Weber 1995) and in the general case the solution is given in terms of the 
Fourier integrals based on Bessel functions with index v = \f— b + (1 — a) 2 /4 that also defines the index of 
equation (50). In the specific case of the dust-dominated FLRW cosmological model the equations have three 
different indices v — 3/2, v — 5/2, and v = 7/2. Fourier integrals constructed from the Bessel functions 
with such indices allow one to find solutions of equations (50) explicitly in terms of the retarded integrals. 
Without going into specific details of our calculations we give the final results: 



1. General solution of equations (50) with index v = 3/2: 



F(r),x) 
as) 



1 d 



ri&n \ v ) 



f d 3 x' p-\x-x'\ 
J \x - x'\ ^ ~ ~ J US ( U ' X '^ du ' 



2. General solution of equations (50) with index v = 5/2: 



F( V ,x) = 



d_ 

drj 



1 d 

i] drj \rj 

d 3 x' . 
Tift 



5/2 



rt]-\X-X'\ rll 

'I) / udu / S{v,x')dv . 

Jrto Jun 



IX — X 



\X — X 



3. General solution of equations (50) with index v = 7/2: 



F{n,x) 

*(r),x) 



d f 
= ri— i 



1 d 



drj \rjdr] 
d 3 x' 



1 d 



rj drj \rj 



\x — x 



77 (V~ \x-x'\) 



+ F 7 / 2 , 

T)-\X-X'\ 



udu 



I vdv I 

Jun Jvq 



S(w,x' 



dw 



Here T v (j\, x) is a general solution of the homogeneous form of equation (50), 
T v {r),x) = Re J f„(k)(kr))<- 1 - a)/2 J v (krj) exp(ik ■ x)d 3 k , 



(51) 
(52) 

(53) 
(54) 

(55) 
(56) 

(57) 



where J v is the Bessel function of order v, k = \k\, and f v (k) is the complex Fourier amplitude of T v {rj, x). 
It is important to realize that the homogeneous solution T v (j], x) could in general involve free gravitational 
waves which would not be completely eliminated by gauge transformations. 
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